Dimer state of spin-1 Bosons in an optical lattice 
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In this paper we consider spin-1 Bosons, such as 23 Na, trapped in an optical lattice, in the regime 
of one particle per site. We argue that the ground state is expected to be the dimer phase in one, 
, two, or three dimensions, thus realizing a state that has so far been studied only theoretically. 

PACS numbers: 03.75.-b, 03.75.Lm, 03.75.Mn 
Tremendous progress has been made recently on trapping and cooling atoms. Greiner et al [1] has succeeded in 
observing the Bosonic Mott-Hubbard transition. In the experiment, the Bosons involved are 8 Rb atoms and the 
lattice potential is provided by the standing waves of three orthogonal laser beams. For weak lattice potentials, the 
^ 1 system remains a superfluid (Bose-Einstein condensate). When the lattice potential is increased beyond a critical 
strength, the tunneling rate for atoms between the different wells becomes weak compared with the repulsion for two 
■ atoms residing in the same well, the system enters into the Mott regime. In this case the number of atoms in each 
well is essentially fixed at an integer value (here, one). Phase coherence and hence superfluidity is lost. [2,3] 

The 87 Rb atoms in the experiment of Greiner et al [1] were hyperfine-spin-polarized. In this paper, we discuss the 
Qh| interesting physics that can be realized if 23 Na atoms are employed instead of 87 Rb and if the atoms are not polarized. 
Bose-Einstein condensation of unpolarized 23 Na has already been achieved by the MIT group [4]. The 23 Na atoms 
have hypcrfine spin (hereafter simply as "spin") 1 in their lower energy manifold, and the interaction among them 
is antiferromagnetic [5,4] (in contrast to 87 Rb, where it is ferromagnetic). We shall confine ourselves to the Mott 
regime where there is essentially one atom per cell, and discuss the arrangement of the hyperfine spin states (or, more 
precisely, the projections) for the atoms in this lattice. We shall argue that the ground state is expected to be the 
dimer phase in one, two, or three dimensions. This is in contrast with spin-1 electronic systems where the effective 
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^ ■ Hamiltonian is the antiferromagnetic Heisenberg model. There in one dimension, the system would be in the Haldane 
J>^\ phase, whereas long range Neel order is expected to be present in two and three dimensions. [6] 

Consider then a (cubic) lattice formed by three orthogonal laser beams, with effective potential for the 23 Na atoms of 
— '! the form V(x, y, z) = V^sin 2 (fca:) + VySivP^ky) + V z sm 2 (kz) where k = 2£ is the wave- vector of the lasers. The strength 
Vx,y,z of the sinusoidal potentials are proportional to the intensities of the laser beams and can be adjusted separately. 
As in Ref [3] we assume that only one orbital state is involved for each well. The effective Hamiltonian for our system 
OO can be written in a Bose-Hubbard form [3] (generalized to Bosons with spin). The hopping matrix elements, t x ,y,z 
are dependent on the directions of hopping. The order of magnitude of t x is given by \t x \ ~ (ErVx) 1 / 2 &~ 2 ( Vx ^ Eb ^ 1/2 
where Er = h 2 k 2 /2m and m is the mass of the atoms. Though we are considering one particle per well, the system 
can exist in excited states where the wells are multiplely occupied. We shall confine ourselves to the regime where 
' these energies are large compared with the hopping energies |t x ,j/,z|- Thus we need only consider those excited states 
with at most two particles per well. The extra energy is described by the Hubbard repulsion Us which depends on 
the total spin S of the two particles involved. Us ~ {EnYxVyVzY^kas where 05 is the scattering length in total spin 
£ = 0,2 channel. Excited states with S = 1 are not allowed due to the identity of the Bosons with one orbital state 
C/5 ' per well. 

We next construct, in the standard manner, the effective Hamiltonian for the spins in the subspace of exactly one 
particle per well, assuming \t\ « t/0.2 (<< the excited energies ~ {Esy x ,y,zY^ of higher orbital states in the wells: 
this latter inequality is typically satisfied since as « A [3]). Let us first consider two sites, labelled by 1 and 2. For 
t — > the energy is independent of the spin configurations. For finite but small t, we can perform perturbation in t. 
It is convenient to classify the states according to the total spin Stot of the two sites. To second order in t, the energy 
is lowered due to hopping by the amount — 4t 2 /Ua,2 for total spin Stot = 0, 2 [7] but zero for total spin 1. With Si,2 
the spin-1 operators for sites 1, 2 and using Si • S2 = — 2, — 1, 1 for Stot = 0, 1, 2 respectively, we find that the effective 
Hamiltonian can be written as 

Hi2 = e + J(S X • S 2 ) + K(Si ■ S 2 ) 2 (1) 

where J = -f^, K = - |^ and e = J - K . For 23 Na atoms, U 2 > U > [5,4,8], hence K < J < 0. As a 

matter of fact, since a 2 ~ 52a^ and a « 46a^ (here = Bohr radius), (£/ 2 — Uq)/Uq = (a 2 — a )/a << 1. Hence 
\J\ ~ \K\ and I J - K\ « \K\ or \ J\. [For 87 Rb, U > U 2 > 0, in that case J < K < ]. 

Thus the Hamiltonian (1) differs significantly from the Heisenberg one [9] familiar in ionic crystals. There, although 
the total spin of an ion can be larger than 1/2 and hence (Si • S 2 ) and (Si ■ S 2 ) 2 are independent operators, the 
Hamiltonian can, to a good approximation (ignoring spin-orbit interactions etc), be written simply as Hnei = J(Si ■ 
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S2). For perturbation to second order in the hopping of electrons between the two ions, there cannot be terms in 
Si • S2 with powers higher than 1 in the effective Hamiltonian, since two hoppings of electrons with spin 1/2 can only 
change the z component of the spin of an ion by at most 1, whereas, e.g., (Si • S2) 2 consists of terms which can change 
that component by ±2. Higher order perturbation in hopping can give rise to terms higher powers in (Si • S2) but 
only with much smaller coefficients when hopping is small compared with the Hubbard repulsion. 

Another peculiar fact about Hamiltonian (1) is also apparant. If the spins were classical vectors, K < J < would 
require that the two spin vectors be parrallel in the lowest energy state. However, since U2 > Uq > 0, Es=o < ^S=i 
and hence the spins on the two sites actually prefer to be anii-ferromagnetically correlated. The strong quantum 
mechanical nature of the spins will be of significance below. 

Let us now begin with the case where V x « V VtZ , so that one can ignore couplings along y and z directions. We 
thus have a collection of one-dimensional spin-1 chains. Ignoring the e term in eq (1) not of relevance below, our 
effective Hamiltonian is thus (for one chain) 

H = Y, [ J(Si-S l+1 ) + X(S ; -S i+ i) 2 ] (2) 

= \J J 2 + K* [ cos ^ ( S * • S '+i) + sin ^ ( S * • S '+i) 2 ] ( 3 ) 

where the sum is over the site labels I and the second relation defines 7. For K < J < — < 1 < — §• If 
\U 2 - U \ « U a <U 2 as for 23 Na, 7 -> -^f , whereas if U 2 » U > 0, 7 -» -§. [For J < K < as in 8 ^Rb, 

-7T< 7 <-f]. 

Hamiltonian of the form in (3) has been studied before, with most efforts in the region where the ground state is 
expected to be in the Haldane phase ( — j < 7 < j). [10] There has been much fewer studies in the range of 7 of 
relevance here. Though it is agreed that for 7 near — ^ the ground state should be in a dimer state [11], results for 
other 7's have been somewhat controversial, especially near — . In particular Chubukov [12] suggested that there 
is a critical 7 C ( — ^ < j c < — |) such that the dimer state is unstable for 7 < j c , where instead the nematic phase 
should exist. A later study [13] claims otherwise. Recently Demlcr and Zhou [14] considered possible ground states 
for spin-1 Bosons in optical lattice as in the present paper, did not include the dimer state in their discussion. 

We would like to first re-address this issue by variational ansatz. We shall denote the three possible spin projection 
states at a given site by |+ >,|0 > and |— >, and write the wavefunctions by specifying these states at each site. The 
nematic phase is given by 

V Nem = I ... 0000 ... > (4) 

This state is in direct analogy with the corresponding "polar" phase [5] in the bulk. This state has < S; >= and 
< S?, x >=< Sf y >t^< Sf z >. [15] For the dimer state, the usually employed ansatz is [10,12] 

* Dimer = ■■■ *12 *34 •■• (5) 

where 

*i2 = 4=(| + ->+!-+> -|00 >)l2 (6) 
v3 

is a singlet ( 5tot = 0) pair formed by sites 1 and 2 (the subscripts label the sites). This state is shown schematically in 
Fig 1 (a). (There is another state degenerate with (5) with all pairs shifted by one lattice site. ) In the non-interacting 
spin-wave approximation, the the S = 2 modes become unstable at j c = — 7r + tan _1 | w — 0.667T. Chubukov [12] then 
concludes that for — ^f- < 7 < j c < — |, (with j c possibly renormalized) , the dimer state is unstable, and further 
speculates that the correct ground state should be the nematic state. 

Here we shall revisit this question of stability by studying an improved variational ansatz. We shall write again eq 
(5) but with the pair wavefunction given by 

*i2 = --/== (| + ->+|- + >-C|00>) 12 (7) 

etc., with ( our variational parameter. If ( = 1, our ansatz reduces to that in eq. (6). Notice that if \(\ — > 00, then 
^Dimer ~ * ^Nem- It is simple to evaluate the expectation values of H. We find, with ( = |C|e 4< ^, 
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< #12 >= j^tt^ {-(1 + 2|C|cos0)cos 7 + (3 + 2|C|W + |C| 2 )sin 7 } (8) 



■ICI 

for the bond (the part of H) between sites 1 and 2, and 



< &3 >= [2+ | C|2]2 {(3 + 2|C| 2 + |C| 4 )sin 7 } (9) 
for the bond between sites 2 and 3. The energy E per site is thus given by 

E Dl mer = ^ [< #12 > + < #23 >] (10) 

The cj) dependent part of E arises only from < H\2 >, and is proportional to (sin 7 — cos 7 )cos</> with a positive 
coefficient. For our region of 7, (sin 7 — cos 7 ) < 0, hence the energy is minimized at <fi = 0, i.e., £ real and positive. 
Henceforth we shall put <j> = and restrict < ( < 00. 
By expanding eq (10) near ( = 1, we find 
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E D irner = |^-C0S7 + ^Sin 7 J + ^ ^COS 7 - ^Sm 7 J (( - if + ... (11) 

Thus £ = 1 is a relative energy minimum provided 7 > 7 C = — n + tan _1 | w — 0.667T defined above. This result is 
in accordance with the discussions following eq (6). However, we here do not interprete this as an instability of the 
dimer state, but rather that a better or more generalized ansatz is required as is done here. Eoimer as a function of 
Q is plotted in Fig 2 (a). For 7 < j c we find that the minimum energies occur at 1 < ( < 00. Thus though the pair 
wavefunction differs from the singlet state, the system is still dimerized. (It can be seen easily that < #12 >y^< #23 > 
provided ( 7^ 00.) ( increases without limit when 7 decreases towards — ^f. To examine the stability of the dimer 
state versus the nematic phase, we expand eq (10) as a function of 77 = C^ 1 near ( = 00 (77 = 0). We find 

E D imer = 2 SUI7 + 2(shl7 - COS7)t7 + ... (12) 

where the ellipsis means terms of order r/ 2 or higher. Note that the energy of the nematic phase is given by putting 
77 = 0, i.e., E^em = 2sin7. Thus the nematic phase is never stable for our 7 region of interest, since (sin7 — COS7) < 0. 
The behavior of E Dimer as a function of £ for large £ is plotted in Fig 2 (b). 

Since (sin7 — COS7) changes sign at 7 = — , the nematic phase (4) is more stable than the dimer state (7) for 
7 < — However, for — 2n < 7 < — the ferromagnetic state 

^Ferro = | ... + + + ... > (13) 

[or any other state obtained by applying the lowering operator S toti _ arbitary (limited only by twice the number of 
sites) number of times], with energy per site Ef erro = COS7 + SU17, becomes more stable than the nematic phase. Thus 
the nematic phase has no regime of stability for the Hamiltonian in eq (3). 

The ansatz wavefunction (7) does not have definite total spin S to t unless £ = 1. Moreover, the wavefunction should 
be improved by including spin-wave fluctuations. However, we expect that the projection onto definite Stot [16] and 
correcting the states with correlations between pairs will not change qualitatively the picture given above. 

Consider now finite V y but still V z — ► 00. We thus now have a collection of two dimensional spin-1 lattices. In 
general the tunneling matrix elements along the x and y directions, hence the strength of the effective spin interactions, 
can be unequal. The effective spin Hamiltonian for our lattice is of the same form as eq (2), except that we need 
two labels (l x ,l y ) for each the lattice point, and the interactions can be between nearest neighbor pairs along both 
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the x and y directions, with interaction constants J x , y , K x y . With t y = 6 y t x , we have J y = S y J x and K y = S y K x . 
Without loss in generality we can take < 5 y < 1. Since the excitations of the dimer phase are gapped [10,12], we 
expect that the dimer phase is stable at least for small S y . It is straight-forward to generalize our ansatz (5) (7) 
above to the two dimensional case, with schematic wavefunction as shown in Fig 1 (b). For each site, there is one 
"strong" bond along x-dircction with energy as in eq (8), one "weak" bond also along x with energy as in eq (9), and 
two "weak" bonds along y direction with total energy given by 25 y times eq (9). From this total energy, we can see 
that our discussions for one dimension is basically unaffected. In particular since eq (9) does not contribute any term 
linear in 77 = £ _1 , the discussions below eq. (12) is qualitatively unchanged. Within our ansatz, the nematic phase is 
still always unstable towards the dimer phase for any — ^ < 7 < — | and < S y < 1. 
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Under the above variational ansatz, the states are degenerate with respect to the spatial arrangements of the pairs. 
In particular, the ansatz wavefunctions corresponding to Fig 1 (b) and (c), with the pairs forming a rectangular 
and triangular lattice respectively, are completely degenerate in energy. However, this degeneracy will be lifted once 
spin- wave fluctuations are taken into account. For ( = 1 the spin wave spectrum can be found as in ref [12]. There are 

V ( ~*\ " ^" / ^ " V ( ~~*\ " 1 / ^ 

three S = 1 modes and five S = 2 modes, with dispersions LOi(q) = X± 1 and uj 2 {q) = X 2 1 H — 

where X\ = cosj — 3sin7, X 2 = 3(cos7 — sin7), Y\ = |(2cos7 — sin7) and Y 2 = |sin7. X 2 , Y~i, Y 2 are all 

positive in our region of 7). Denoting the distance between sites by a = ir/k, for the rectangular lattice, u Tect (q) = 
cos(2q x a) + 28 y cos(q y a), and for the triangular lattice, v trl (q) = cos(<f • a) + S y cos(q ■ (a — bj) + S y cos(q ■ b). Here 
a = (2a)x and b = a(x + y) are the lattice vectors for our triangular lattice. The correction to the energy from the 
spin- waves is given by AE = X^ll [ w i(<7) — (-^i ~~ ^l^io)] + § [^2(9) — (-^2 + ^2^(5)]}. Evaluating this energy for 
small S y , we find that the rectangular lattice has lower energy. Assuming that the lattice type applies for the entire 
region of 7 of interest here, we conclude that the ground state for two dimensions should be as shown in Fig 1 (b). 
Similar considerations suggest that, for three dimensions, the lattice for ground state should be tetragonal. 

The distinguishing property of the dimer state is the doubling of the unit cell, while the spins do not have long 
range order. Thus the periodicity of the excitations with period Aq = Tr/a = k = 2ir/\ for a certain direction (x 
above) would be the signature of the dimer state. These excitations can in principle be detected by scattering. 

In conclusion, we have pointed out that the trapped Bosonic atoms provide the opportunity to realize a quantum 
mechanical state predicted in theory of quantum magnetism so far not tested experimentally. Obtaining this state in 
the laboratory would further widen our the play-ground for quantum many-body systems 
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FIG. 1. Dimers. Thick lines represent the pairs such as eq (7). 
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